Quasiparticle Thermal Conductivities in a Type-II Superconductor at High Magnetic 

Field 
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We present a calculation of the quasiparticle contribution to the longitudinal thermal conductiv- 
ities Kxx{H,T) (perpendicular to the external field) and Kzz{H,T) (parallel to the external field) as 
well as the transverse (Hall) thermal conductivity hixy{H,T) of an extreme type-II superconductor 
in a high magnetic field [Hci <C -ff < Hci) and at low temperatures. In the limit of frequency and 
temperature approaching zero (il — > 0, T — > 0), both longitudinal and transverse conductivities 
upon entering the superconducting state undergo a reduction from their respective normal state 
values by the factor (F/A)^, which measures the size of the region at the Fermi surface containing 
gapless quasiparticle excitations. We use our theory to numerically compute the longitudinal trans- 
port coefficient in borocarbide and A- 15 superconductors. The agreement with recent experimental 
data on LuNi2B2C is very good. 

PACS numbers: 74.25Fy, 74.60Ec, 74.70Dd, 4.70Ad 



I. INTRODUCTION 



The low-temperature, high- field region in the H-T phase 
diagram of an extreme type-II superconductor is the 
regime where the Landau level quantization of the elec- 
tronic energies within the superconducting state is well 
defined, i.e. the cyclotron energy huc > A, T, F where 
A = A(T, H) is the BCS gap, T is the temperature and 
r = F(a;) is the scattering rate due to disorder. This 
regime should be contrasted with the more familiar op- 
posite limit of low magnetic fields and high temperatures 
where electrons occupy a huge number of Landau levels 
and where the temperature and/or impurity scattering 
broaden these levels and reduce the significance of Lan- 
dau quantization. With the Landau level structure fully 
accounted for one discovers a qualitatively new nature of 
quasiparticle excitations at high fields: for fields H be- 
low but near Hc2 the quasiparticle spectrum is gapless 
at a discrete set of points on the Fermi surface. These 
gapless excitations reflect a coherent quasiparticle prop- 
agation over many unit cells of a closelyi-jpacked vortex 
lattice with fully overlapping vortex coresoHa. The pres- 
ence of such low-lying excitations makes an s-wave, "con- 
ventional" superconductor in a high magnetic field some- 
what similar to an anisotropic, "unconventional" super- 
conductor with nodes in the gap. In the low-temperature 
and high-field regime however, the nodes in the gap re- 
flect the center- off-mass motion of the Cooper pairs in 
the magnetic field, in contrast to d-wave superconduct- 
ing cuprates where such nodes are due to the relative 
orbital motion. This gapless behavior in bulk systems 
is found to persist to surprisingly low magnetic fields 
H* ~ 0.2 — 0.5Hc2- Below H* the gaps start opening 
up in the quasiparticle spectrum and the system eventu- 
ally reaches the low-field regime of localized states in the 
cores of isolated, well-separated vorticesH. At the present 



time the strongest evidence for the quantization of quasi- 
particle orbits within the superconducting state comes 
from observations of the de Haas-van Alphen CdHvA) 
oscillations in various superconducting material^. The 
persistence of the dHvA signal deep within the mixed 
state, with the frequency of oscillations still maintaining 
the normal state value, can be attributed to the pres- 
ence of a small portion of the Fermi surface containing 
gapless quasiparticle excitations, surrounded by regions 
where the gap is largeOa. 

Another useful probe of low-energy excitations in su- 
perconductors is measurement of their thermal transport. 
The simultaneous measurements of the field dependent 
longitudinal Kxx{H, T) and transverse KxyiH, T) thermal 
conductivities are now feasible experimentally and can 
yield information on both quasiparticle dynamics and the 
pairing mechanism. The dependence of transport coeffi- 
cients on magnetic field is currently a hotly debated issue 
in the scientific community in light of the experimental 
observation of field independent plateaus in the longitu- 
dinal thermal conductivity of high temperature |^uper- 
conductors (HTS) at low fields (iJ^i < < i?c2)Q. The 
field-independent attributed to the d^^^y^ pair- 

ing mechanism at low fields and to tlieiiodal structure of 
the resulting quasiparticle excitationsoEi. The presence of 
propagating gapless quasiparticles in the superconduct- 
ing state at low temperatures and high magnetic fields 
should also lead to transport properties qualitatively dif- 
ferent from those found in s-wave superconductors at low 
fields, where the number of thermally activated quasipar- 
ticles is exponentially small and the only contribution to 
the thermal conduction is found along the field direction 
and originates from the bound states within vortex cores. 
Recently, the thermal conductivity of the borocarbide su- 
perconductor LuNi2Bap was measured down to T = 70 
mK by Boaknin et aZJlJ in a magnetic field perpendicular 
to the heat current from iJ = to above Hc2 = 7 Tesla. 
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In the limit of T — > 0, a considerable thermal transport 
was observed in the mixed state of the superconductor 
{Hci < H < Hc2), indicating the presence of delocal- 
ized low-energy excitations at the Fermi surface. On the 
other hand, no thermal transport was observed at zero 
field, a result consistent with the s-wave superconducting 
gap without nodes at the Fermi surface. 

The purpose of this work is to examine the contribu- 
tion of low-energy quasiparticles to the thermal transport 
of conventional, i.e. extreme type-II superconductors in 
the regime of high magnetic fields and low temperatures. 
The paper is organized as follows: In Sec. II we develop 
the Kubo formalism for the transport coefficients within 
the Landau level pairing mechanism while in Sec. Ill we 
incorporate disorder in the Green's function description 
of a three-dimensional superconductor in a high magnetic 
field. We use the formalism of Sec. II and III in Sec. IV 
to examine both longitudinal Kxx{H, T) and Kzz{H, T) as 
well as transverse KxyiH,T) conductivities in the regime 
of the freqency ^ and the temperature T — + 0. Fi- 
nally, in Sec. V we report on numerical calculations of 
the thermal transport in the borocarbide LuNi2B2C and 
the A-15 superconductor VsSi and compare our theoret- 
ical plots with the available experimental data. 



II. KUBO FORMALISM IN THE LANDAU 
LEVEL PAIRING SCHEME 

Thermal conductivities can be calculated within the 
framework of the Kubo formalism as |a-|linear response of 
a system to an external perturbationliil 
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where Ulf{n) = Uij{in n + iS) and 



(1) 



U,j{in) = J d^xid^X2ll,j{l,2;in), 

% (1,2; if]) = - / dre^f^- < T,j,(xi, r)j,(a;2, 0) > (2) 
Jo 

is the spatially averaged, finite temperature thermal 
current-current correlation function tensor and f3 = 
l/ksT. In order to derive the heat current operators 
j(l) and j(2) at the space-time point 1 = (xi,r) and 
2 = (x2,|Q^ we follow the standard s-wave derivation in 
zero fieldEj and generalize it to our case of a non-uniform 



gap at high fields. A similar approach was recently uti- 
lized in Ref. 9 for a d-wave superconductor in zero field. 
The heat current carried by the quasiparticles can be 
computed within the standard variational procedure as 



dC 



dC 



from the Lagrangian density 



(3) 



:(V4VV^„-VV'iVa)-A 



-^SV'lV'laV'-aV'a (4) 

where all the energies are measured with respect to the 
chemical potential and ip = dip/dt. The Nambu's two- 
component field operators '0a = are written in a 
compact notation for the sake of brevity. In extreme 
type-II superconductors, as soon as the magnetic field 
satisfies H ^ Hci{T) the vector potential A = A(r) can 
be safely assumed to be entirely due to the external field 
H = V X A. This holds over most of the H ~T phase 
diagram. We have used a simple BCS-model point in- 
teraction V{tx — V2) = —gb(v\ — v-i) in expression (||). 
To lowest order in the concentration of impurities, the 
electron-impurity interaction can be omitted in comput- 
ing the heat current. The effect of disorder will be in- 
cluded later in the Green's functions for a superconduc- 
tor. The variational procedure (0) yields 



1 

~2to 
eA 



d 



(V 



eA' 



■7rT(V+ . 

OT^ CI 



(5) 



for the heat current operator. With this definition 
it is straightforward to calculate the correlation tensor 
n.y (1, 2; zfi) within the usual Hartree-Fock approxima- 
tion (i.e., bare bubble approximation) defined by 

< T,^,(l)Vfe(2)^|(2')V'l(l') >^ G„(l,2')Gfe,(2,l') 

(6) 

where G^i is the Nambu's matrix Green's function. In- 
serting (||) into (||) and with the help of (^) the current- 
current correlator becomes 



-i(f7 + /i)(Vi 



eA(l') 



eA(l) 



eA(2), . ,^ eA(2'); 
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Cl 



Tr [t3G(1, 2', n + /i)T3G(2, 1', 2^2' 



(7) 
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where rs is a Pauli matrix, fi = 27rm//3 are bosonic Mat- 
subara frequencies and 1 = xi . 

It was shown in Ref. 1 that the mean-field Hamiltonian 
corresponding to the Lagrangian density (^) can be diag- 
onahzed in terms of the basis functions of the Magnetic 



Sublattice Representation (MSR), characterized by the 
quasi-momentum q perpendicular to the direction of the 
magnetic field. The eigenf unctions of this representation 
in the Landau gauge A = y, 0,0) and belonging to 
the m-th Landau level, are 



^fc^ ,q,m 



(r) = 



exp [i{q, 



^A_exp(*fc.C)E^-P(*^fc' 
irk , 



ikqyby) X 



— )x - l/2{y/l + qj + —im„.{^ + (fe + 



a 



(8) 



where ^ is the spatial coordinate and is the momen- 
tum along the field direction, a = (a, 0) and b = {b^, by) 
are the unit vectors of the triangular vortex lattice, 
I — yhc/eH is the magnetic length and LxLyLz is the 
volume of the system. Hm{x) is the Hermite polyno- 
mial of order m. Quasimomenta q are restricted to the 
first Magnetic Brillouin Zone (MBZ) spanned by vectors 
Qi = iby/P,-bJP) and Q2 = (0,2a/Z2). 

Normal and anomalous Green's functions for a clean 
superconductor in this representation can be constructed 



G2i(l,2;u;)= C^fc,,-q(l)Cfe.,q(2)^n ^)(9) 



where w = (2m -I- 1)tt/(3 are the electron Matsubara fre- 
quencies. Similar expressions can be written for the re- 
maining two Nambu matrix elements. In writing (j^) we 
have taken into account only diagonal (in Landau level 
index n) contributions to the Green's functions. This 
is a good approximation in high magnetic fields where 
A/huJc ^ 1 and the number of occupied Landau levels 
He is not too large, which is the case for the extreme 
type-II systems under consideration. In this situation 
we are justified in using the diagonal approximatioroB, 
in which the BCS pairs are formed by electrons belong- 
ing to mutually degenerate Landau levels located at the 
Fermi surface while the contribution from Landau levels 
separated by hujc or more is included in the renormaliga- 
tion of the effective coupling constant (g g{H,T))t3. 
As long as the magnetic field is larger than some critical 
field H* (T) the off-diagonal pairing does not change the 
qualitative behavior of the superconductor in a magnetic 
field. The critical field H* at T ~ can be estimated 
from the dHvA experiments to be ~ Q.BHcs, for A-15 and 
~ 0.2Hc2 for borocarbide superconductorsQ. 

When Nambu matrix (^ is inserted in Eq. ^ and 
the space average in (^) is performed, the longitudinal 
Hxxii^) = ^yyi'i-^) a-nd transverse (Hall) Ilxy{ifl) = 



-Ilyx{iil,) current-current correlation functions become 



71+1 



Tr[T3Gnikz,q, + ?;cj)r3G„+i (fc^, q, iu) 
±T3Gn+i{kz,(l,in + iuj)T3Gn{kz,(l,iuj)] (10) 

where the -I- sign corresponds to U^xii^), the — sign cor- 
responds to illxy{i^) and oj = {2m+l)Tr/ jS are electronic 
Matsubara frequencies. On the other hand the longitu- 
dinal (parallel to the external magnetic field) current- 
current correlation function Hzz{i^) becomes 

"^^^^"^ = 4^^^ ^ fcf(f^ + 2^)^ 
xTr[T3Gnikz,(l,in + iuj)T3Gnikz,(l,iuj)] . (11) 

In order to perform the summation over the Matsubara 
frequencies uj, we introduce a spectral representation for 
the Nambu matrix G„(fc2,q, w) as 



/oo 
-OC 



An{kz,q,uji) 



1— (12) 

tUJ — UJl 



where the spectral function matrix An{kz, q, to) is defined 
as 



A„(fc^,q,cj) 



1 



/mG;^*(fc„q,c^). 



(13) 



When the spectral representation of the Green's func- 
tions (^ is used back in (^ and (^ respectively, we 
obtain 



n+l 



xTr[T3An{kz,q,uJi)T3An+i{kz,q,uj2) 

±T3An+likz,q,UJl)T3Anikz,q,LU2)] X 5 (14) 



and 



n,kz ,q 



xTr[T3An{kz,q,uji)T3An{kz,q,uj2)] x S (15) 



4 



where S contains Matsubara sums i.e. 



{ifl + iuj — uji){iuj — UJ2) 



(16) 



The sum can be evaluated in the standard way by pick- 
ing up the., contributions from each of the poles of the 
summancO. After the analytic continuation ifl ^ fl + iS 
we obtain the retarded function Sret 



072 ~ ^1 + + «<5 

where np{uj) is the Fermi function. 



(17) 



In order to obtain the imaginary part of 11^ (fl) we 
need to find an imaginary part of 5*^6* when calculating 
the longitudinal conductivity KxxiS^^T) and Kzz{^iT). 
On the other hand, since Sret enters the expression for 
illxyi^) in ( |l^ ) we need to find the real part of —Sret- 
Using the identity 



= P inS{x) 



X + i6 X 



(18) 



and taking the imaginary part of ([l^), we find that the 
diagonal conductivities become 



l^xx ^yy 



+00 



duo- — :; — — — ^ X 



2^2 



n + 1 



Tr [T3A„(fc2, q, uj + f7)T3A„+i (fc^, q, u) + T3yl„+i(fc^, q, uj + rj)r3An(k2;, q, u)] 



(19) 



and 



n,kz ,q 



+00 



— rti?([j + SI) 



y2 



Tr [T3A„(fcz,q, w + 12)r3^„(fcz,q,w)] 



Similarly, taking the real part of (^^ with the help of 
(O) yields for the off-diagonal conductivity: 



J 



(20) 



Kxy Kyx 1 [^°° ^ {2uj + npi^o) — npiuo + Vl) 



T2 







n + 1 



Tr [r3S„+i(A:z,q,w + 0)r3A„(/cz,q,cj) - T3Bn{kz,q,uj + f^)T3yl„+i(fcz, q, t 



(21) 



where the function Bn{kz, q, cj + fi) is defined as 



III. GREEN'S FUNCTIONS IN THE PRESENCE 
OF DISORDER 



/+00 
dull 
-oc 



An{kz,q,uji) 

UJ — UJl + fl 



(22) 



Before further discussing expressions (|l^), (|20| ) and 
( ^ ) we should go back to the question of spectral func- 
tions or, alternatively, Green's functions for the super- 
conductor in a magnetic field. The Green's function for 
the clean superconductor can be easily found following 
Ref. 14 with their "Fourier transforms" in the quasi- 
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momentum space expressed in the Nambu formahsm as 



I f iLU + e„{kz) -A„„(q) 

(zc^)2 - E„{k,,q) \ -A,*„(q) - e„(fc,) 



where 



(23) 



En,p{kz,q) = phujc ± + I A„+p,„_p(q)|2 



2m 



+ hujc{n + 1/2) - fi (24) 



is the quasiparticlc excitation spectrum of the supercon- 
ductor in high magnetic field near points k^ — ±kpn — 



2m{fj. — hiOcin + 1/2)) /h^ calculated within the diago- 
nal approximatioiiBB, where A/huc ^ 1- For the quasi- 
particles near the Fermi surface {k^ ~ kpn) it suffices to 
consider only the i?n,p=o bands. The gap, A„„(q), which 
in the MSR representation can be written as 



A„m(q) 



(-1)" 



%/2 2^+"Wn\m\ 



Eexp(in—k^) x 
n 



nk 



nk , 



exp{2ikqyby - (g^; + — )^^^)i?„+„jV2(gfa; + (25) 

turns to zero on the Fermi surface at the set of points in 
the MBZ with a strong linear dispersion in q. The exci- 
tations from other, p 7^ in (|2^), bands are gapped by 
at least a cyclotron energy and their contribution to the 
quasiparticle transport can be neglected at low tempera- 
tures (T ^ A(T, H) <^ hujc)- Once the off-diagonal pair- 
ing in (^ is included, the excitation spectrum cannot be 
written in the simple form ( p^ ) and a closed analytic ex- 
pression for the superconducting Green's function cannot 
be found. Nevertheless, when these off-diagonal terms 
are treated pertubatively as in Ref. 3 the qualitative be- 
havior of the quasiparticle excitations, characterized by 
the nodes in the MBZ, remains the same. This state- 
ment is correct in all orders of the perturbation theory 
and therefore is exact as long as the pertubative expan- 
sion itself is well defined i.e. as long as H > H*{T). 
Once the magnetic field is lowered below H* , gaps start 
opening up at the Fermi surface signalling the crossover 
to the low-field regime of quasiparticle^tates localized in 
the cores of widely separated vorticesLJ. 

In a dirty but homogenous superconductor with the 
coherence length ^ much longer than the effective dis- 
tance ^imp over which the impurity potential changes 
(f/Cimp >> 1), the superconducting order parameter is 
not affected by the impurities and still forms a perfect 
vortex lattice. For such a system, the bare Green's func- 
tion in (H) is dressed via scattering through the diagonal 
(normal) self-energy E^(iix')_-and off-diagonal (anoma- 
lous) self-energy i;^„(q, icj) JlJ A dressed Green's func- 
tion is obtained by replacing lo with a) and A„„(q) with 
A„„(q) in (|) where 



iuj = iLU — (iuj) 
A„„(q) = A„„(q) + E;^„(q, lu). 



(26) 



In order to calculate the spectral functions in ( |l9| ) 
the analytical continuation should be performed so 
that G'ret(fcz, q, t^) = G{kz,q,iuj lo + iS) where 



(CO) 



-^N.A 



{iuj 



LO + i6) with the impurity 



scattering rate in the superconducting state defined as 
r{u!) — —ImS^^iiijo). It was shown by the authors in 
Ref. 14 that the anomalous self energy docs not quali- 
tatively change the form of the gap function A„„(q) at 
low energies and therefore E;^„ (q, uj) will be neglected in 
further calculations. At the same time, the real part of 
the normal self energy Yj^ (uo) can be either neglected or 
absorbed into Cnikz). 



IV. THERMAL CONDUCTIVITIES IN T ^ 
LIMIT 

We are interested in calculating thermal conductivities 
in (^ and (|2^) in the limit of and small T such 
that T < T{lj). In the limit of f2 ^ the difference of 
Fermi functions in (n9[) becomes 



npii^ + ^) — nF{to) dnp 



(27) 



This function is sharply peaked around w = at very 
low temperatures so that we are justified in expanding 
the integrand in (|l^) and ( ^l]) around a; = up to 
second order in to and setting the scattering rate to a 
constant F = r(ct; = 0). In the high- field superconduc- 
tors the largest contribution to the thermal conductivity 
comes from the quasiparticle excitations at the Fermi sur- 
face with momenta q such that A(q) < maa:(T, F) while 
the excitations gapped by large A(q) give exponentially 
small contributions. Therefore, in order to simplify the 
integration over the MBZ and summation over the Lan- 
dau level index in ( [l9| ) and (|2^) we linearize the excitation 
spectrum (|2^) around nodes at the Fermi surfaca3. This 
enables us to obtain approximate but analytic expresions 
for thermal transport coefhcients which capture the qual- 
itative behavior near Hc2. Keeping thisin mind and with 
the help of the Sommerfeld expansionlia the longitudinal 
conductivity Kxx{H,T) — Kyy{H,T) and Kzz{H,T) at 
low temperatures become 



and 



KzziH,T) _ fvy , 7^2 /^^y 



K^,iH^O,T) yAj 15 V A 



(28) 



(29) 



K^^{H,T) is the thermal conductivity of the normal 
metal in a magnetic fieldEJ 



{H,T) 



3 2TO*r (hiOc)^ + 4r2 



T 



(30) 
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where = ^Fn- '^^ electronic density in the sys- 

tem. On the other hand the transverse conductivity 

K 



xy\ 



{H,T) = ~Ky^{H,T) becomes 



"1^ 



5 TT 



Kxy{H,T) 

(31) 

where K^y{H,T) is the off-diagonal thfirmal conductivity 
of the normal metal in magnetic fielcii3 



(rr ^ (^^^)' r 



(32) 



Relations (gSD, (|29|) and (|3|) obtained when T/A < 1 
within the "linearized spectrum aproximation" tell us 
that there is still a considerable thermal transport in the 
mixed state of the superconductor. This is in stark con- 
trast to the exponential suppresion of transport charac- 
teristic of an s-wave superconductor in zero field. Fur- 
thermore, relations (p8|), ( p9| ) and ( |3l] ) indicate that when 
passing from the normal to the superconducting state, 
both longitudinal and transverse transport coefficients 
k/T are reduced from their respective normal state val- 
ues by the factor ^ (F/A)^ (the term linear in (T/A)^ is 
negligible at low temparatures even for very clean super- 
conductors). The factor ~ (F/A)^ measures the fraction 
of the Fermi surface Q containing gapless quasiparticle 
excitations at T = 0. The size of Q is determined by 
both the total number of nodes in the excitation spec- 
trum ( p^ and the areas in different branches where the 
BCS gap A is very small but not necessarily zero. This 
result, obtained here for the thermal coefficients, is con- 
sistent with the behavior of some other superconducting 
observables that measure the presence of low-energy ex- 
citations at the Fermi surface. One such experimentally 
confirmed behavior is the reduction of the de Haas-van 
Alphen (dHvA) oscillation's amplitude in both A-15 and 
borocarbide supierconductors when the sample becomes 
superconducting3. The drop in the overall amplitude in 
passing from the normal to the superconducting state re- 
fiects the presence of a small portion of the Fermi surface 
~ Q containing coherent gapless excitations while the rest 
is gapped by large A.eI 



V. COMPARISON WITH EXPERIMENT 



Recently, the longitudinal thermal conductivity of the 
borocarbide superconductor LuNi2B2C-|Was measured 
down to T = 70 mK by Boaknin et aZJ13 in a magnetic 
field from _ff = to above Hc2 — 7 Tesla. In the limit of 
T — > 0, a considerable thermal transport is observed in 
the mixed state of the superconductor [Hd < H < i?c2), 
indicating the presence of delocalized low-energy excita- 
tions at the Fermi surface. The authors argue that this 
observation is strong evidence for a highly anisotropic 
gap function in LuNi2B2C, possibly with nodes. On the 
other hand, no sizable thermal conductivity was observed 
in zero field, the result expected for a superconducting 
gap without nodes. 

The Boaknin et al. result is consistent with the obser- 
vation of dHvA oscillations down to fields H* ~ i?c2/5 
in YNi2B2C, a close cousin of LuNi2B2C as well as 
in VsSi where, the oscillations persist down to fields 
H* ~ i7c2/2.u It was shown by us that the drop in 
the dHvA amplitude observed in these experiments can 
be atributed to the quantization of quasiparticle orbits 
within the superconducting state which results in the 
formation of nodes in the gapcl. This quantum regime 
behavior in fields H* < H < Hc2 is due to the center- 
off-mass motion of the Cooper pairs, in contrast to the 
d-wave or anisotropic s-wave where nodes in the gap are 
due to the relative orbital motion. Therefore it makes 
sense to compare the theory developed in this paper 
with the experimental data in Ref. 10. A quick check 
tehs us that the expression (|8|), where T{H)/A{H) < 1, 
does not hold through the entire range of fields used in 
the experiment. Therefore, we numerically compute the 
longitudinal thermal conductivity directly from Eq. 19, 
without using any additional aproximations, for both the 
borocarbide superconductor LuNi2B2C as well as for the 
A-15 superconductor VsSi. 



In the limit of fl 
yields 



and T ^ the expression ( p^ ) 



J 



12771 TlLUr 



(" + i)E 



(33) 



r 



where the number of Landau levels involved in super- 
conducting pairing tt-c = Ep/Tiujc varies as a function 
of magnetic field. In the borocarbide superconductor 
LuNi2B2C Tie can be estimated as ~ 33 at i?c2 a-nd 
77c ~ 1147 at field oi H = 0.2 Tesla (these numbers 



were obtained using an effective mass of 0.3577ie and 
Fermi velocity vp — 2.76 x 10^ cm/s as reported in 
Ref. 18). On the other hand, the number of occupied 
Landau levels in the A-15 superconductor VaSi is much 
larger: Uc ~ 241 at Hc2 = 18.5 Tesla and Uc ^ 4470 
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at _ff = 1 Tesla (we used an effective mass of 1.7me 
and Fermi velocity vp — 2.8 x 10^ cm/s from Ref. 19 
in this estimate). The scattering rate F = r(w = 0) 
in (^) is, in general, modified relative to the normal 
state scattering rate Fq when the system becomes su- 
perconducting. Indeed, the self-consistent calculation of 
F in Ref. 14 gives T{H) = ^FoA(i7)/2. We assume 
A{H) = A^l - H/Hc2 which is a good approximation 
for the range of fields used in the experiment. 

The dashed line in Fig. 1 shows the magnetic field de- 
pendence of the quasiparticle thermal conductivity k/T 
for the borocarbide superconductor LuNi2B2C in the 
limit of T ^ obtained by the numerical evaluation 
of (^11), where values for the BCS gap A and normal 
state inverse scattering rate Fg are taken from Ref. 10. 
The full circles-|in Fig. 1 are the experimental data of 
Boaknin et a/.tJ for the same superconductor. The full 
line in Fig. 1 shows the theoretical plot obtained by nu- 
merical evaluation of ( p3| ) for the A- 15 superconductor 
VsSi where values for A and Fq are taken from Ref. 19. 
There is a significant difference in the behavior of k/T 
in these two superconducting systems characterized by 
much smaller thermal transport in VaSi when compared 
to the transport in LuNi2B2C in the same range of mag- 
netic fields. This observation indicates that the number 
of gapless or near gapless excitations at the Fermi surface 
in the VaSi is very small at H Hc2- In order to under- 
stand this difference, one has to notice that in magnetic 
fields H < H* = 0.5i?c2 the number of occupied Landau 
levels in VsSi system is huge (< 4500) and that VaSi is 
away from the regime of coherent gapless excitations of 
high fields. On the other hand, the number of occupied 
Landau levels in LuNi2B2C is much smaller (< 1000) 
and it seems that there are still many gapless excitations 
left at low fields. It is suprising though that significant 
thermal transport exists down to if ~ 0.015iJc2, a field 
much smaller than the critical field H* = 0.2Hc2 for this 
system, where most of the quasiparticle spectrum should 
be gapped. Note, however, that a possible source for 
such a significant transport at low fields might be the 
highly anisotropic s-wave gap function in LuNi2B2C, as 
suggested in Ref. 10. If there is such anisotropy along one 
or more directions on the Fermi surface, the range of va- 
lidity of our theory may bepSxtended to fields lower than 
the simple estimate for H* .Ej In this regard, we alert the 
reader that at the lowest fields in Fig. 1 {H ^ 0.04iic2) 
our theory is stretched to its very limits and its quanti- 
tative accuracy diminishes. 



VI. CONCLUSIONS 

In this paper we develop expressions for the longitu- 
dinal and transverse quasiparticle thermal conductivities 
for an extreme type-II superconductor in a magnetic field. 
We utilize the Landau level formalism of superconduct- 
ing pairing in a magnetic field to obtain, within the Kubo 
mechanism of linear response to an external perturba- 
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FIG. 1: Magnetic field dependence of the quasiparticle lon- 
gitudinal thermal conductivity computed from Equation (^) 
for LuNi2B2C (dashed line) and VaSi (full line^ Full circles 
represent experimental data of Boaknin et alE3 The verti- 
cal dotted line indicates the normal-superconducting transi- 
tion at Hc2 = 7 Tesla for LuNi2B2C. The upper critical field 
for VsSi at Hc2 = 18.5 Tesla is not shown. For LuNi2B2C 
we have used experimentally determined values for A = 4.4 
meV and Fq — 0.5A from Ref. 10 as well as the effective 
mass m* = 0.35me from Ref. 18. For VsSi, A = 2.6 meV, 
Fo = 0.61 meV and m* — l.Ym^ were taken from Ref. 19. 



tion, thermal currents perpendicular and parallel to the 
external magnetic field. From there, current-current cor- 
relation functions are introduced within the Matsubara 
finite temperature mechanism in order to derive closed 
expressions for thermal conductivities Ky (r2,T). We ex- 
amine the transport coefficients Hij/T in the limit of 
il ^ and T ^ and find that there is considerable 
thermal transport in the mixed state of a superconductor 
with an s-wave symmetry due to the creation of gapless 
exitations in the magnetic field. This is in contrast to 
the zero field thermal transport which is exponentially 
small for an s-wave superconductor with no nodes in the 
gap. Furthermore, when passing from the normal to the 
superconducting state the thermal coefficients Kij/T be- 
come reduced with respect to their normal state values 
by a factor ~ (F/A)^ which measures the fraction of 
the Fermi surface that contains coherent gapless or near 
gapless excitations in a magnetic field. In this respect, 
thermal conductivities behave similarly to the dHvA os- 
cillations in which the amplitude is also reduced at the 
superconducting transition. Finally, we numerically com- 
pute the longitudinal thermal conductivity for two realis- 
tic superconducting systems, the borocarbide LuNi2B2C 
and the A-15 superconductor VaSi. The thermal trans- 
port in LuNi2B2C is much larger in magnitude than the 
thermal transport of VaSi at the same field. This re- 
sult indicates that the borocarbide LuNi2B2C might still 
be in the regime of delocalized quasiparticle states even 
at fields much lower than the critical field H* 0.2Hc2 
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(estimated from the dHvA experimentsHuj) . The agree- 
ment of our theoretical plot with the-experimental data 
for LuNi2B2C taken by Boaknin et alS^ over a wide range 
of fields used in the experiment is suprisingly good. 
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